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Abstract 

This work calculates the first correction to the leading order qq dilepton production rates due 
to shear viscosity in an expanding gas. The modified rates are integrated over the space-time 
history of a viscous hydrodynamic simulation of RHIC collisions. The net result is a hardening 
of q± spectrum with the magnitude of the correction increasing with invariant mass. We argue 
that a thermal description is reliable for invariant masses less than M max ~ (2tqTq) / (rj / s) . For 
reasonable values of the shear viscosity and thermalization time M max ~ 4.5 GeV. Finally, the 
early emission from a viscous medium is compared to emission from a longitudinally free streaming 
plasma. Qualitative differences in </j_ spectrum are seen which could be used to extract information 
on the thermalization time, viscosity to entropy ratio and possibly the thermalization mechanism 
in heavy- ion collisions. 
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I. INTRODUCTION 



There is a general consensus that the early matter produced at RHIC behaves as a near 
perfect fluid [l|. This conclusion was born out of the success of ideal hydrodynamic descrip- 
tions 2, [J of both hadron transverse momentum spectra and elliptic flow measurements 
up to 1.5-2 GeV/c. Although it is too early to draw any definitive conclusions most likely 
the deviations from ideal hydrodynamic behavior can be ascribed to dissipative effects. 
This has already been suggested in some of the recent works on dissipative hydrodynamics 

B,y,y,H,y,y,[ioi. 

In addition to hadronic observables which interact strongly and therefore depend only 
on the final state of the medium, electromagnetic probes are emitted throughout the entire 
space-time evolution reaching the detector without any final state interactions. In terms of 
computing observables there is a big difference, since the resulting transverse momentum and 
elliptic flow spectra depend only on the final freezeout hypersurface, whereas the resulting 
dilepton yields depend on the full space-time volume. A consistent description of heavy-ion 
phenomenology should use the same space-time evolution for both hadronic spectra and 
dilepton observables. 

In this work we calculate the first viscous correction to dilepton emission from quark 
anti-quark annihilation in a dissipative medium. The kinematic region when a thermal 
description is reliable is found by requiring that the viscous corrections are small. When 
the viscous corrections become large a kinetic description is really required. The viscous 
rates are then integrated over the space-time history of a hydrodynamic simulation of RHIC 
collisions. We show how shear viscosity modifies the transverse momentum and invariant 
mass spectrum. We find that the inverse slope of the transverse mass spectrum is sensitive 
to both the thermalization time as well as the shear viscosity and can therefore be used in 
order to learn about the early stages of heavy-ion collisions. Finally a comparison is made 
with dileptons produced from a free-streaming quark-gluon plasma. 

II. DILEPTON RATES 

The rate of dilepton emission from a quark-gluon plasma due to qq annihilation is given 
in the Born approximation as 
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dN f d 3 k! d 3 k 2 

where q = (go, q) is the virtual photon's four momentum and M 2 = [E\ + E2) 2 — (ki + k^) 2 
is the photon's invariant mass. Throughout this work we consider massless quarks; therefore 



E\ } 2 = y ki,2 + m g m {k-i^l- The function f(E,T) is the quark or anti-quark momentum 
distribution function, which in thermal equilibrium is given by f(E,T) = 1/(1 + e E ^ T ). 

In the above expression t> 12 is the relative velocity of a quark anti-quark pair and a(M 2 ) 
is the qq cross section. Both expressions are well known from the literature [ll| and are 
given by t> 12 = 2 e*e 2 anc ^ a (,M 2 ) = — ■ The integral over the quarks' momentum can 

be done analytically with the result 



dN a 
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where n± = l/(e^o±lql)/2r + y and f b ^T) = l/(e q °/ T - 1). 

III. VISCOUS CORRECTION TO THE DILEPTON RATES 

In order to account for dissipative effects in the dilepton emission rate we include the first 
viscous correction to the quark and anti-quark's distribution function in eq. [TJ This approach 
neglects any space-time inhomogeneities and assumes that the distribution functions relax 
to their dissipative forms much quicker than the medium evolves. Ideally, one could solve the 
Baym-Kadanoff equations which would take non-equilibrium effects into account. We note 
that the leading order born qq rates do not contain pinch singularities which suggests that at 



least to leading order one may neglect space-time inhomogeneities [12| . This approximation 
allows us to calculate the dilepton emission rates locally in a space-time volume d 4 x. 

As shown in HQ , [ij]] viscosity modifies the ideal distribution function. The resulting 



correction for fermions is 

/(*) ^ /(*) + 2(6 ^ )T2 /( fc )[ 1 - f(k)]k a k^ aP , (3) 

where 7i a/3 = r](V a Up) = r](V a U/3 + Vpu a — |A aj gV p M p ) and r] is the shear viscosity not be 
confused with the space-time rapidity rj s . The coefficient C\ can be computed analytically 
for a massless fermion gas and is given by C\ = 147r 4 /1350C(5) w 0.97. Substituting the 
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above result into the born annihilation rate (eq. [T]) and keeping terms up to first order in 
rj/s (quadratic in momentum) one obtains: 
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In simplifying the above result we have used the fact that the permutation of k\ <-» &2 
gives the same contribution after integration. We write the final result as the sum of the 
ideal and viscous correction 



f- q = h{q) + V^ l2 " {q) ^^ 
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Since I 2 is a second rank tensor depending only on u a and q a it can be decomposed as 



T a/3 



a g a P + a lU a u p + a 2 q a q p + a 3 («V + u^q a ) 



The final result will contain the term I 2 TT a p. By making use of the identities u a Ti a p 



(?) 



9 — only the term with coefficient a 2 will be non-vanishing. a 2 is found by using the 
identity a 2 = P a pl® where the projection operator in the local rest frame of the medium is 



1 
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We now quote the final result for the first viscous correction to the born dilepton annihi 
lation rates: 



4 



dN 
~¥q 



N c a 2 e 2 q 



12tt 4 

where we have defined 



/t(?0 , r)[1 + ^ M ^ )] __|L_ 



(9) 



6 2 (?o,|q|) = t4 r + /(^i,T)/(g -^i)(l-/( J Ei)) 

iqi 



(3g 2 - |q| 2 )£ 2 - 3g ^iM 2 + -M 4 



dE l 
(10) 

and _E± = |(go ± |q|)- For large invariant masses (M/T ^> 1) one can replace the Fermi 
distribution with the classical Maxwell-Boltzmann distribution. In the viscous correction 
to the distribution function this amounts to substituting //(l — //) — > /mb- In this limit 
an analytic expression can be found for the viscous correction to the dilepton rates. In the 
limit that (u ■ q)/T 3> 1 the resulting expression is given as 



dN = N c a 2 el 
d 4 q 12vr 4 
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3(e + p)T^ * " ap \ ' (H) 
where as before C\ ~ 0.97. We find that the above result holds at the accuracy of a few 

percent for M > 3 GeV at T = 400 MeV. 

A feature of the viscous correction is that it does not modify the invariant mass spectrum. 
This is seen by looking at either of the above forms of the viscous correction (eq. [TT]or [9]). In 
going from d 4 q to dM 2 the resulting integral will be a second rank tensor depending on u a 
only. The most general form the result can take is a linear combination of terms proportional 
to g al3 and m"^, which both vanish when contracted with n a/ 3. 

Before performing the full space-time evolution we illustrate the effect of the viscous 
correction using a simple model for the gradients. We consider a ID Bjorken expansion 
without transverse flow. The viscous component of the stress-energy tensor can be easily 
computed and is given as 
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By substituting the above result into eq. [TT] an analytic expression can be found for the 
dilepton yields in the limit that M/T 3> 1. After performing the integration over r) s the 
result is 
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FIG. 1: Kinematic regions where the viscous correction is less then order one. More precisely, the 
boundary is set by the condition \5f /fo\ < 0.8. 

where K u (x) is the modified Bessel function evaluated at x = m±/T. Fig. [T] shows 
the kinematic regions where the viscous correction is small (i.e. dN V i S / 'dJS \deai < 0.8) and 
therefore dictates when using a thermal description of dilepton production is suitable. The 
criterion that dictates when hydrodynamics is applicable can be written as 



and can therefore be separated into a condition on the medium, r]/s, and a condition 
on the experimental setup, 1/(tT). Throughout this work we always set f]/s = 0.2. The 
region surrounded by the solid line is for 1/(tT) « 2.2 corresponding to a temperature of 
450 MeV at r = 0.2 fm/c. The region surrounded by the dotted line is for l/(rT) m 0.65 
corresponding to a temperature of 300 MeV at r = 1 fm/c. At earlier times the viscous 
correction is larger and the allowed region is smaller. 

The results shown in fig. [T] should only be taken qualitatively. Transverse flow alleviates 
the situation, opening up the boundaries shown above. Also, the result presented was in 
the limit M/T » 1 where an analytical result was obtained. Fig. [TJ is still useful, since 
it still serves as a qualitative picture where the viscous corrections become large even after 
including flow and relaxing the classical limit. Outside of the kinematic boundaries a thermal 
description may no longer be reliable. 

Throughout the remainder of this work we now resort to eq. M, which is evaluated nu- 
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FIG. 2: (color online) Dilepton transverse momentum spectra for T = 0.4 GeV and n/s = 0.2 at 
r = 1 fm/c for a boost invariant expansion with no transverse flow. 

merically, in order to compute the dilepton yields accurately at all masses. Fig. [2] shows the 
dilepton spectrum generated for a temperature T=0.4 GeV at proper time r = 1 fm/c and 
using a viscosity to entropy ratio of rj/s = 0.2. 

Figures of the invariant mass spectrum are not shown because, as discussed earlier, the 
spectrum is unmodified when including the viscous correction. Looking at the transverse 
momentum s pec trum a hardening of dileptons is seen that is reminiscent of the single particle 
spectrum in [131 ] . The magnitude of the viscous correction is dictated by both r//s as well 
as the proper time. At earlier times the shear between the longitudinal and transverse 
directions is larger resulting in bigger corrections at smaller proper times due to the 1/r 
factor in eq. 



IV. EVOLUTION MODEL 

In order to model the space time evolution of the collision region we use the results of jj], 
which is summarized in this section. The hydrodynamic model is a 1+1 dimensional boost 
invariant expansion with initial conditions tuned in order to simulate Au-Au collisions at 
RHIC energies (\/s = 200 GeV). Dissipative corrections to the ideal hydrodynamic expansion 
is treated using a second order relaxation scheme first proposed by Ottinger and Grmela 

Detailed plots of the hydrodynamic solution with and without viscosity is shown in (J). 
For modest values of the shear viscosity (r]/s < 0.3) dissipative effects did not integrate to 
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give large changes to the ideal hydrodynamic solution. The net effect of a finite viscosity 
was twofold. First, the longitudinal pressure is initially reduced causing a slower decrease 
of energy density per unit rapidity at early times. The reduction of longitudinal pressure is 
accompanied by a larger transverse pressure which drives larger transverse velocities. The 
larger velocities then cause the energy density to deplete faster at later times. 

Even though the changes to the ideal hydrodynamic result is small a full viscous simu- 
lation is still needed in order to have access to the velocity gradients which enter into the 
dissipative corrections of the quark and anti-quark distribution functions. 

The hydrodynamic model is started at tq = 0.2 fm/c in order to account for some of the 
pre-equilibrium production of dileptons which will contribute at larger masses. Dileptons are 
produced as long as the temperature of the medium is greater than a critical temperature 
taken as T c = 0.170 GeV. We do not look at dileptons produced during a mixed phase or 
hadronic phase in this work. At any space-time point the hydrodynamic model provides 
the three independent terms of the stress tensor (ii rr ,r 2 7i^ and T 2 7i vv ). The equation for 
g a g /3 7r Q , / 3 used in eq. M is given as 

q a qP 7 r ct p = g 2 cos 2 (0)tt"" + q\ sin 2 (0)rV^ + mi sinh 2 (^)r 2 7r w 

+mi cosh 2 (r] s )v 2 iT rr — 2m± cosh(r] s )q± cos(9)vK rr . (15) 

Fig. [3] shows the resulting transverse momentum spectrum after the full space-time 
integration at two invariant mass points: M = 0.525 GeV (left) and M = 2.625 GeV 
(right). First the red curve shows spectra generated from an ideal hydrodynamic simulation 
(r]/s = 0). Next the green curve shows the spectra generated from a viscous simulation 
having rj/s = 0.2 but without including the viscous correction to the distribution function. 
This curve therefore shows the effect that viscosity has on modifying the ideal hydrodynamic 
equation of motions. We find that a finite viscosity leads to a slight increase in the overall 
yield. This is due to the fact that a finite viscosity causes the energy density to deplete more 
slowly at early times. This effect therefore brings about an effective increase in the lifetime 
of the simulation above the critical temperature. We find ~ 30% increase in the low mass 
region and ~ 50% increase in the higher mass region. 

Finally, the blue curves in fig. [3] show the viscous result including the viscous correction 
to the distribution function. We find that the magnitude of the viscous correction increases 
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FIG. 3: (color online) Dilepton transverse momentum spectra after the full space-time integration 
of a boost invariant expansion with arbitrary transverse expansion and azimuthal symmetry. 

with the invariant mass. This was similarly observed in fig. [1] where the range in having 
viscous corrections of order less than one (as shown by shaded regions) decreased in size with 
increasing mass. The simulation results are discussed in more detail in the next section. 



V. DISCUSSION 

In order to further understand the viscous corrections, the effective temperature (T e ff) 
of the dilepton spectrum from the full space-time integration is constructed as a function of 
invariant mass. The effective temperature is found by fitting the transverse mass spectrum 
at a given mass to the following expression, 

'"' V oce^/^. (16) 



dM 2 m±dm±dy 

In this work the fit is done in the transverse momentum region of 0.5 < q±(GeV) < 2.0. 
As expected, we find that the transverse mass spectra does not exactly fit the above form. 
Actually, other ranges in could have been chosen where the fit works better. However, 
the results are qualitatively the same and therefore a different choice in q±_ range will not 
change the discussion that follows. If a quantitative comparison were to be made with data, 
it would be more appropriate to compare to the actual q± spectra instead. Regardless, T e ff 
still serves as a useful quantity since it probes the average temperature of the medium as 
well as the radial flow profile and viscous correction. Looking at fig. [3] we expect the viscous 
correction to increase the effective temperature, with larger corrections at higher masses. 
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FIG. 4: (color online) Effective temperature as a function of invariant mass. 

In fig. 0] the effective temperature is shown as a function of invariant mass. The solid blue 
curve labeled ideal shows T e ff for an ideal (j]/s = 0) hydrodynamic expansion started at 
proper time r = 0.2 fm/c. The shape of the curve is dictated by the underlying radial flow 
as well as the average temperature of the emission region. At higher masses the dominate 
source of dileptons is from the higher temperature regions, which occur at earlier proper 
times. This explains the slight rise in T e ff with mass. We now look at the solid red curve, 
which is generated from a hydrodynamic evolution having r]/s = 0.2. In this case we do not 
include the viscous correction to the distribution function and the resulting modifications to 
the effective temperature are due to changes in the hydrodynamic evolution. As discussed 
earlier, modifications from viscosity to the hydrodynamic solution are small and we therefore 
don't expect to see large deviations from the ideal case. This is indeed the case. 

We now focus the discussion on the role of the viscous correction to the distribution 
function. The green curve in fig. H] shows the effective temperature of dileptons coming 
from a viscous medium having 77/s = 0.2 from a simulation started at a proper time of 
t = 0.2 fm/c. The result is that the effective temperature increases greatly as a function 
of invariant mass. From the magnitude of the correction, the upper bound of the domain 
of hydrodynamics is found to be at most M max « 2.0 GeV for this parameter set. There 
are two reasons why the viscous correction increases with mass. First, there is the explicit 
mass dependence in the viscous correction itself. This is easiest to see by looking at the 
approximate form, eq. [E2 The second reason is because the high mass contribution is mainly 
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produced in the early, high temperature stages of the evolution. Looking at eq. [T2]the viscous 
correction grows like 1/r at early times. In order to see the effect of the early emission a 
final simulation is done (dashed-black curve) where the hydrodynamic evolution is started 
at To = 1 fm/c. In this case the viscous corrections are more modest and M max « 4.5 GeV. 

It is therefore a combination of both the thermalization time as well as the magnitude of 
r]/s that dictates when a hydrodynamic description is reliable. Since the effective tempera- 
ture rises so quickly with mass, as long as there is non- vanishing viscosity, there will always 
be a mass region outside of the region of a hydrodynamic description. From eq. [13] one can 
extract an approximate condition for the mass. Since most of the particle yield is at low q± 
we should guarantee that the viscous correction is small at q± = 0. Furthermore at high and 
intermediate masses the ratio of Bessel functions is approximately one. Then the condition 
that the viscous correction must be less than of order one can be expressed as 

M max » ^ . (17) 
rj/s 

When the viscous corrections to the spectra become large a kinetic approach is re- 
quired. One can ask whether the viscous correction at early times mock up the effects 
of off-equilibrium production that would be taken into account by a full kinetic theory. In 
order to test this hypothesis T e ff spectra is calculated from a free streaming non- interacting 



gas of quarks [171 Il8| | . We should point out that our treatment is very similar to the re- 



cent work of [19j. In this model the initial parton distribution is taken as thermal with 
the temperature chosen in order to reproduce the thermal dilepton number given by the 
hydrodynamic simulation. Starting with the thermal initial condition at r = 0.2 fm/c the 
total dilepton yield is found by integrating the free streaming result IA10I up to a final time 
of r = 1.0 fm/c. The details of this calculation is given in the appendix. We now discuss 
the results. 

We consider two scenarios. The first is running the hydrodynamic simulation starting 
at ro = 0.2 fm/c until T c . The second scenario runs the free streaming model from 0.2 < 
r(fm/c) < 1.0. Then at 1.0 fm/c the hydrodynamic evolution is started and ran until T c . 
We should stress that the second model is not very realistic since the free streaming model 
is not asymptotic with the hydrodynamic evolution at r = 1 fm/c. A future work might use 
a more realistic model for the evolution of the distribution function then the proposed free 
streaming case. For example, one could start with an initially anisotropic distribution which 
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FIG. 5: (color online) Left: Dilepton invariant mass spectra. Right: Effective temperature as a 
function of invariant mass. 



evolves to its thermal form from multi-quark scattering 20||, at which point a hydrodynamic 
evolution is started. However, we expect the true result to lie between the two scenarios 
used in this work. 

Figure [5] shows both the invariant mass spectrum (left) and effective temperature (right) 
for the two scenarios outlined above. The curves to compare are the hydrodynamic simu- 
lation started at To = 0.2 (labeled Hy. r = 0.2 fm/c) and the sum of the hydrodynamic 
simulation started at To = 1 and the free-streaming (labeled FS+Hy.). We first note that 
the invariant mass spectrum is qualitatively the same for the two scenarios and it would 
not be possible to discern between the two scenarios from experimental data. Qualitative 
differences do appear in the T e ff spectrum. First we find that the free-streaming with hydro 
solution mimics the early time hydro only solution at low masses. However, at high masses 
the two result diverge when the viscous correction can no longer be trusted. While the 
early hydrodynamic solution increases greatly with mass the free streaming solution flattens 
off at higher mass. We therefore argue that through a detailed analysis of q± spectra one 
could hopefully extract information on the thermalization time, viscosity to entropy ratio 
and thermalization mechanism in heavy-ion collisions. 



VI. CONCLUSIONS 



In conclusion, we have calculated the first viscous correction to dilepton production from 
leading order qq annihilation. The rates are then integrated over the space-time history of 
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a viscous hydrodynamic simulation of RHIC collisions. We argue that a thermal description 
is only reliable for invariant masses less than w (2roT 2 ) / (77 / s) and above this a kinetic 
description is required. For rj/s = 0.2 and r = 1 fm/c this corresponds to M < 4.5 
GeV. We have shown that viscosity does not change the invariant mass distribution but 
strongly modifies the transverse momentum distribution and can therefore be used to extract 
information on both the viscosity to entropy ratio as well as the thermalization time. Finally, 
we have also made comparisons with an initially free streaming QGP. Qualitative differences 
in transverse momentum are seen, which could again possibly be used to learn about the 
thermalization mechanism. 
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APPENDIX A: DILEPTON YIELDS FROM A FREE STREAMING QGP 

In this appendix the dilepton yields are derived for a free streaming boost invariant 
expansion. The starting point is the collision-less Boltzmann equation 

P^f(p,x) = 0, (Al) 

where f(p,x) will be considered as the phase-space distribution for the quark and anti- 
quark. Under the assumption of boost invariance as well as homogeneity in the transverse 
plane the Boltzmann equation can be written as 

3 T /-& x / = 0, (A2) 

where x = V ~ Vs- The initial condition is such that the quark distribution is isotropic 
and starts from local thermal equilibrium, f(p, r = r ) = eTO /r +1 - One can write po = u-p = 
p± cosh(x) by using the assumption of boost invariance and homogeneity in the transverse 
plane. We note that even for quarks out of equilibrium it still might be useful to use the 
equilibrium form of the distribution function where T is instead considered as an effective 
temperature describing the initial state. The solution of eq. IA2I at any time r is 
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f(p, x) = . (A3) 

e ^V 1+sinh2w (^) + i 

With the explicit form of the distribution function available one can calculate the dilepton 
rates using the same kinetic theory expression used before (see eq. [I]) 

= I I ^^ / ^ 1 ' x)/(p2 ' x)t;i2a(5(4) ^ 1+P2_g) - (A4) 
Making use of the expressions for the relative velocity and cross section as quoted earlier 
the above equation can be expressed in the following form 

dN B j d'x J d i p 1 d%6(p 2 1 )6(pl)6^( Pl +p 2 - q)f(pi,x)f(p 2 ,x) , (A5) 



d 4 q 



327ra 2 e 2 



where B - -^p-- 
First let us quote some well known identities: 



d 4 x = rdrdr] s d 2 x± = 7iR 2 Tdrdr] s 
d 4 p = -dp 2 dy p p±dpj_d(p p 
d 4 q = MdMq±dq±dyd(j) 
6®{P -q) = 45(P 2 - M 2 )6(y p - y)6(P L - g ± )5(0 P - 0) 

where P M = P\ + p 2 and y and </> are rapidity and angle in the transverse plane. We 
place a subscript P on quantities to indicate they are derived from P M . The free streaming 
dilepton rate can now be expressed as 

(IN f f 

— = Bjd 4 'xj dy i p lj ±dpx t ±d<p 1 dy 2 p2,±dp2,±d^ 2 f{pi,x)f(p2,x)x 

5(P 2 - M 2 )5(P 2 ± - ql)5(y P - y)5(<p p - 0) (A6) 

Since the distribution function is boost invariant the integral over rj s is trivial due to the 
delta function, 5(yp — y). By defining y± = y x ± y 2 and <f>± = 4>i ± 02 the delta functions 
can be rewritten as 
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After rewriting the integration variables as dy\dy 2 = dy+ 2 dy and 
integral over 0_ and y_ can be done explicitly yielding 
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The delta function in the above equation enforces the following constraints 
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2pi,_Lf>2,± 

Let us make a further shift of variables, p± = p\ t ± ±p2,±- The constraints ( 1A9I) then take 
particularly simple form. The final expression is 



dN 



dM 2 dydq 2 



N a 2 e 2 

487T 4 



V/(M 2 + qj - pl)(M* + qi- pi) - qjjtfl - p 



dy+ I dp + dp_(p 2 + -pi)x 

Qx J -tlx 

1 /(pi,r)/(pa,T) (A10) 



where 
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f(pi,r) 



1 + exp 



/(p 2 ,r) 



1 + exp 



y_ = sinh 



-i 



2y/(M 2 + qj - pj) (M 2 + qj £g) 



(AH) 
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